tspan = [3.9 , 4];
y0 = [8, 2];
[t,x] = ode45('odefunc', tspan, y0);


